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ABSTRACT 

We performed ^^C0(1 - 0), ^'^00(1 - 0), and HCN(1 - 0) interferometric observa- 
tions of the central region (about 450 pc in radius) of M82 with the Nobeyama Mihimeter 



O . 

. Array, and have successfully imaged a molecular superbubble and spurs. The center of 

the superbubble is clearly shifted from the nucleus by 140 pc. This position is close 
Qh| to that of the massive black hole (BH) of > 460 Mq and the 2.2/zm secondary peak 

O I (a luminous supergiant dominated cluster), which strongly suggests that these objects 

' may be related to the formation of the superbubble. Consideration of star formation in 

■ the cluster based on the infrared data indicates that (1) energy release from supernovae 

can account for the kinetic energy of the superbubble, (2) the total mass of stellar-mass 
^ . BHs available for building-up the massive BH may be much higher than 460 Mq, and 

^ , (3) it is possible to form the middle-mass BH of 10^ — 10^ Mq within the timescale of 

the superbubble. We suggest that the massive BH was produced and is growing in the 
intense starburst region. 

Subject headings: black hole physics, galaxies: individual (M82, NGC 3034), galaxies: 
ISM, galaxies: starburst, ISM: bubbles 
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1. INTRODUCTION 

A starburst is a star formation event with a high star formation efficiency. This usuaUy 
occurs in the galactic nuclear /central regions. Starbursts and successive supernova explosions 
cause drastic changes in kinematic and physical conditions of interstellar medium, and produce 
(expanding) superbubbles, chimneys, and/or large-scale outflows. There is observational evidence 
for these structures in hot and/or atomic gas in many galaxies (e.g. Fabbiano 1989; Tcnorio-Tagle & 
Bodenheimcr 1988), but there is much less clear evidence in molecular gas. Starburst phenomena in 
the nuclear regions are also thought to be closely related to active galactic nuclei (AGNs) powered 
by massive black holes (MBHs; more massive than a stellar-mass BH), such as fueling AGNs and/or 
causing the growth of MBHs. One question, "can starbursts produce the seeds of MBHs or middle- 
mass BHs?," is still largely unresolved. 

Nearby (3.25 Mpc; Sandage & Tammann 1975) irregular galaxy M82 (NGC 3034) has kilo- 
parsec scale bipolar outflows which can be seen by optical emission lines (e.g. Shopbell & Bland- 
Hawthorn 1998) and X-ray (e.g. Bregman et al. 1995) (Fig. la). This outflow is believed to be 
made by frequent supernova explosions at the central region in consequence of starburst. Recent 
X-ray studies have found evidences of an MBH with a mass of > 46OM0 in the starburst region 
(Matsumoto k, Tsuru 1999; Ptak & Griffiths 1999). In this region, however, strong dust absorption 
prevents detailed optical/near-infrared observations, so we performed millimeter-wave observations 
which would not be affected by the dust. 



2. OBSERVATIONS 

ApcrtTirc synthesis observations of the central region of M82 were carried out in the "'^■^CO, 
^^CO, and HON J = 1 ^ lines (rest frequency = 115.271 GHz, 110.201 GHz, and 88.632 GHz, 
respectively) with the Nobeyama Millimeter Array (NMA) during 1997 November - 1999 March. 
All the images were obtained using three configurations of six 10m antennas which are equipped 
with tunerless SIS receivers (Sunada et al. 1994), and the system noise temperatures in the single 
side band were about 950 K, 850 K, and 450 K at 115 GHz, 110 GHz, and 89 GHz, respectively. As a 
back-end, we used an XF-type spectro-correlator Ultra Wide Band Correlator (UWBC; Okumura 
et al. 2000), with a total bandwidth of 512 MHz over 256 channels for the -"^^CO observations 
(corresponding to 1300 km s^^ bandwidth with 5.2 km s^^ velocity resolution), and with a total 
bandwidth of 1024 MHz over 128 channels for the ^'^CO and HCN observations (corresponding to 
2800 km s~^ and 3500 km s~^ bandwidths with 22 km and 27 km s~^ velocity resolutions, 
respectively). The band-pass calibration was done with 3C273, and 0923-1-392 was observed every 
10 minutes as a phase and amplitude calibrator. The flux scale of 0923-1-392 was determined by 
comparisons with Mars and Uranus, and has an uncertainty of ~ 20%. 
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3. EXPANDING MOLECULAR SUPERBUBBLE 

The overall distribution in the ^^CO total integrated intensity map shows diffuse spurs north 
and south (minor axis direction) of the galaxy, in addition to the previously identified (Shen k. 
Lo 1995) three prominent peaks (the so-called North-East [NE] and South- West [SW] lobes and a 
central peak). These spurs trace the dark filaments which can be seen in optical B-band images 
(Alton et al. 1999) (Fig. lb). Similar structures have been detected further out in the filaments 
(Nakai ct al. 1987; Kuno &: Matsuo 1997; Alton et al. 1999) which are connected to the large-scale 
Ha and X-ray outflow (Shopbell &: Bland-Hawthorn 1998; Bregman et al. 1995) (Fig. la), but ours 
is the first detection close to the disk. 

A position-velocity (PV) diagram of ^^CO data cut along the major axis of the galaxy (Fig. 2b) 
shows an arc-like deviation from rigid rotation between the central peak and the SW lobe. Such 
deviations are common features of expanding H I superbubbles (Deul k. den Hartog 1990). In 
addition, recent literature also suggests that this deviation implies the existence of a molecular 
superbubble (Neininger et al. 1998; Wei/3 et al. 1999; Wills et al. 1999). However, any shell struc- 
tures have not been detected in spatial image so far, and Wci/3 et al. (1999) concluded that the 
superbubble had already broken toward the minor axis direction of this galaxy. We made a channel 
map, binning over the velocity range Vlsr = 118 to 212 km s~^, which corresponds to the range 
deviates from the rigid rotation (Fig. Ic). The map clearly shows a shell-like structure, with a 
diameter of ~ 210 pc x 140 pc 14" x 9") elongated toward north-south direction, between the 
two peaks. This structure is also visible in other molecular lines such as ^^CO (Fig. Id) and HCN 
(Fig. le). These images indicate that the superbubble is still not broken out of the galactic disk. 
Around the superbubble, however, there are some spurs connecting with the shell structure, so that 
it may be possible that some parts of the superbubble were already broken. 

If we assume that the shell structure is the result of edge-brightening (Wills et al. 1999), the 
southern part of the shell should be moving perpendicular to the line-of-sight with nearly the rigid 
rotation velocity, which corresponds to the systemic velocity of the molecular superbubble. Indeed, 
the velocity field of the southern shell is very close to rigid rotation (Fig. 2c). In addition, the 
PV diagram shows that almost all of the molecular gas except that of the superbubble is on the 
rigid rotation (Fig. 2b). These reasons lead us to conclude that the expansion velocity of the 
superbubble is the largest velocity deviation from the rigid rotation seen in the PV diagram; the 
resultant expansion velocity is ~ 100 km s"^. However, we cannot reject the possibility of an 
expansion velocity of ^ 50 km s"-*^, which is the mean velocity deviation in the PV diagram (Wei/3 
et al. 1999), so we use the velocity range of 50 — 100 km s~^ for the following calculations. Using 
the velocity and the size of the superbubble derived above, the elapsed time from the explosion can 
be estimated as ~ (1 — 2) x 10^ years. Wei/3 et al. (1999) also derived similar timescale of 1 x 10^ 
years. 

Next, we will estimate the energetics of the superbubble. Using the intensity ratios between 
the ^^CO line and the ^^CO and HCN lines, we calculated the CO-to-H2 conversion factor based 
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on the Large- Velocity-Gradient (LVG) approximation (Sakamoto 1999). The resultant conversion 
factor is (1.4ib0.6) x lO^'^ cm~^ (K km s~^)~^, which is consistent with the previous observations of 
~ (1.0-1.2) X 10^0 cm-2 (K km s^^)"! (Wild ct al. 1992; Smith et al. 1991). Hence, we adopted our 
estimated value, and the calculated molecular superbubble mass would be ~ 1.8 x 10^ M©. From 
this molecular gas mass and the expansion velocity, we derive the kinetic energy of the superbubble 
to be ~ (0.5 — 2) X 10^^ erg, equivalent to the total energy of ^ 10^ — 10^ supernovae, and an order 
of magnitude or more larger than that observed in H I superbubbles in other galaxies (< 10^^ erg; 
Tenorio-Tagle &: Bodenheimer 1988). The energy estimated by Wei/3 et al. (1999) was ~ 2 x lO^'* 
erg, similar to ours, although their estimation have a large ambiguity in the volume of the material 
swept up by the explosion. 



4. MASSIVE BLACK HOLE INSIDE THE SUPERBUBBLE 

At the center of the superbubble, there is a hard X-ray variable point source and a 2.2 /xm 
secondary peak (Fig. Ic); hard X-ray observations with ASCA indicate that there is a strong point 
source located in the central region of M82. This source shows time variability in its intensity, 
which indicates the existence of a MBH with its mass of ~ 460 — 2 x 10^ (Matsumoto &: Tsuru 
1999; Ptak & Griffiths 1999). Comparing the ROSAT HRI image (Bregman et al. 1995; Stevens et 
al. 1999) with the false ROSAT image made from ASCA data (the image using similar energy-band 
with that of ROSAT), it was found that the position of the hard X-ray variable source is consistent 
with that of the X-ray peak detected with the ROSAT (Matsumoto k, Tsuru 1999). Detailed 
comparison between the location of this point source and our newly obtained ^^CO map, reveals 
that the X-ray point source is located at the center of the superbubble. Recent high resolution 
X-ray observations with Chandra show that there is a variable source inside the superbubble with 
its peak luminosity of ~ 9 x lO^'^ erg s^^, assuming this source has a same spectrum as the ASCA 
variable source (Matsumoto et al. 2000). This result also supports the conclusion that there is a 
MBH inside the superbubble. On the other hand, the emission from the 2.2 /xm secondary peak 
(Dietz et al. 1986; Lester et al. 1990) seems to be dominated by luminous supergiants (Joy et al. 
1987), which suggest that it is a late phase dense star burst cluster. 

Since these objects are located close to the center of the superbubble, it is natural to think 

that these objects arc related to the superbubble's formation. We therefore discuss the possibility 
that the starburst at the 2.2 ^um secondary peak produces the superbubble and the MBH. We 
first calculated the stellar population of this 2.2 ^m secondary peak cluster assuming an initial 
mass function (IMF), and estimated the starburst evolution as follows: The luminosity of the 2.2 
/xm secondary peak with an extent (full width at half maximum) of 4" (McLeod et al. 1993) is 
equivalent to ~ 1500 M2 supergiants. We assume that stars with initial mass larger than 30 M0 
(short-lived stars with lifetimes < 2 x 10^ yr; Lang 1998) have already exploded as supernovae, 
and that the remaining stars with initial masses of 25 — 30 Mq are now M2 supergiants. Using 
an extended Millar-Scalo IMF (Kennicutt 1983) of dN/dm oc m~^'^ with lower and upper mass 
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limits of 1 and 100, and assuming that there are 1500 stars whose masses are 25 — 30 Mq, the total 
stellar numbers and the total mass formed at the 2.2 /im secondary peak cluster would be about 
8 X 10^ stars and 2 x 10^ Mq, respectively. The number of > 30 Mq stars which would have already 
exploded in this cluster is ~ 4 x 10^, consistent with the estimated number of supernovae needed 
to create the superbubble (10^ — 10^). The observational evidence and IMF calculations suggest 
that the expanding molecular superbubble may be a result of localized starbursts which occurred 
around the position of the 2.2 /xm secondary peak. 

We next discuss the possibility of making a MBH at the starburst in M82. If we assume 
that stars with initial mass of > 25 Mq would create stellar-mass BHs of almost the same mass 
(Brown & Bethe 1994), there would be ~ 4 x 10^ BHs with a total mass of ~ 2 x 10^ Mq in 
the 2.2 iim secondary peak cluster. Assuming an isothermal sphere stellar density distribution 
(Lee 1995), about 4 x 10^ Mq BHs would sink into the cluster center by dynamical friction within 
2 X 10^ yr. This mass is well within the range of that of the MBH estimated from the hard X-ray 
observations. There is also a possibility of star-star merger. If we assume a merging probability 
of 0.1% in the cluster, which is a similar probability to the simulations of stellar mergers in star 
clusters (Portegies Zwart et al. 1999), it is possible to make one 700 Mq star, or several stars of 
a few hundred Mq, in the cluster. Explosions of such very high mass stars can produce > 10^^ erg 
or even 10^^ erg energy hypernovae (e.g. Paczyiiski 1998), which can be seen as 7-ray bursts, and 
may possibly have created the superbubble and MBH. If we set the lower limit of the BH mass, 
Mbh, in this cluster using the hard X-ray observations, and the upper limit using the total mass 
of the stars which seem to be already exploded (> 30 Mq), the range of the BH mass would be 
460 — 2 X 10^ Mq. This mass range suggests that this MBH might be a middle-mass BH. Since 
the cluster mass, Mdustcr, is calculated as 2 x 10^ Mq using the IMF, the ratio of Mbh to Mduster 
would be —3.6 < log(MBH/-Mciuster) ^ —1.0. As shown in Figure 3, this range is on the Mbuige-MBH 
relation for galaxies with supermassive BHs (Magorrian et al. 1998). 

Since there are many other clusters in the central region of M82, another question, "why other 
clusters do not have MBHs?," would arise. One possibility is lack of stellar density. If their densities 
are small, the effect of dynamical friction or star-star merging rate may decrease and it may not 
be possible to make massive objects. The low density clusters would also be affected by tidal force 
(e.g. Taniguchi et al. 2000) and tend to be smaller clusters which are not large enough to make 
massive objects. There is another possibility that the MBH was not created at the cluster but 
supplied from satellite galaxy by minor merger (e.g. Taniguchi et al. 2000). In this case, there are 
two possibilities of the agreement between the position of the MBH and that of the superbubble; 
one is just accidentally overlapping each other, and another is that the superbubble is created by 
the MBH. A BH merging with a compact object (neutron star, white dwarf, etc.) can cause a large 
energy release 10^^ erg; Meszaros et al. 1999). If the MBH meets with a cluster on the way 
sinking toward the galactic center, there are possibilities to merge with some compact objects, and 
as a result, energetic explosions would occur and the superbubble would be made. These discussions 
still have large ambiguities because of lack of detailed information, so that further observations with 
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many wavelengths/frequencies would be needed. 

In future, this middle-mass BH may increase its mass. Since massive stars still exist 10^ 
Mq) at the 2.2 /xm secondary peak cluster, it is possible to feed stars and/or stellar-mass BHs 
to the middle-mass BH. Also, the middle-mass BH is still not at the dynamical center of M82, 
and it may sink down to the center with dynamical friction, where mass is strongly concentrated. 
Therefore there are many ways for the middle-mass BH in M82 to increase in mass and therefore 
grow up to the supermassive BH. Our results may give a new explanation to the reason why some 
of the quasars are embedded in interacting (and therefore active star forming) galaxies, and strong 
7-ray bursts and middle-mass BHs are located either at or offset from galactic nuclei. 

We would like to thank J. Makino, Y. Funato, and anonymous referee for helpful discussion 
and comments. We also thank B. Wallace for carefully reading our manuscript. We are grateful to 
the NRO staff for the operation and improvement of NMA. S.M. and H.M. are financially supported 
by the JSPS Postdoctoral Fellowships for Research Abroad. 
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Fig. 1. — Various scale images of M82. In molecular gas images, value for the colorscale is indicated 

on the top of each figure, and the plus mark indicates the position of the galactic nucleus determined 
from the peak of the strongest 2.2^m. source (Lester et al. 1990). (a) The X-ray image (contours; 
Bregman et al. 1995) overlaid on the Ha image (grayscale; Shopbell &: Bland-Hawthorn 1998). (b) 
^^CO integrated intensity map (contours) overlaid on optical B-band image (grayscale; Alton et 
al. 1999). Bright regions on the grayscale are indicated in black and dark regions are in white. 
The contour levels of the ^^CO map arc 5, 10, 15, 20, 25, 30, 40, 50, •• • ,90ct, where I a = 1.64 Jy 
beam~^ km s~^ [= 23 K km s^^]. The synthesized beam (2'.'8 x 2'.'3 or 42 pc x 35 pc) is shown 
at the bottom- left corner, (c) The ^^CO molecular supcrbubble image, binning over the velocity 
range of 118 - 212 km s"^ The contour levels of the ^^CO map are -10, 10, 20, 30, • • • ,110c7, 
where 1 cr = 8.9 mJy beam~^ [= 127 mK]. The synthesized beam size is the same as (b). The 
diamond mark and the filled circle indicate the central positions of the 2.2|Ltm secondary peak 
(Dietz et al. 1986), and the X-ray point source observed with ROSAT (Stevens et al. 1999) which 
corresponds with the ASCA hard X-ray variable source position (Matsumoto & Tsuru 1999). We 
also indicate the recently observed Chandra X-ray variable source (Matsumoto et al. 2000) with a 
circle, which radius corresponds to its position uncertainty, (d) The HCN molecular superbubble 
image, binning over the velocity range of 134 — 215 km s~^. The contour levels of the HCN map 
are —6,-3,3,6,9,12,15(7, where 1 cr = 3.7 mJy beam~^ [= 39 mK]. The synthesized beam (4''1 
X 3''6 or 62 pc x 54 pc) is shown at the bottom-left corner, (e) The ^^CO molecular superbubble 
image, binning over the velocity range of 127 — 214 km s^^. The contour levels of the -"^'^CO map 
are —3, 3, 6, 9, 12, 15, 18(T, where 1 a = 5.4 mJy beam~^ [= 41 mK]. The synthesized beam (3'.'9 x 
3''4 or 59 pc x 51 pc) is shown at the bottom-left corner. 

Fig. 2. — Position- velocity (PV) diagrams of M82. Values for the colorscale are indicated on 
the top of each figure, (a) ^^CO integrated intensity map. Solid lines indicate the sliced regions 
for PV diagrams displayed below. The plus mark is the same as Fig. lb. (b) PV diagram at 
slice A. The solid line indicates rigid rotation velocity. The region of the superbubble (around 
Q;(B1950)=9''51™43®.6) clearly deviates from the rigid rotation velocity, (c) PV diagram at slice 
B. The solid line indicates rigid rotation velocity. Almost all of the gas at this slice is on the rigid 
rotation velocity. 

Fig. 3. — Correlation diagram between the mass of black holes (Mbh) and those of the host bulges 
or clusters (-Mbuige) in galaxies. The bar mark indicates the estimated mass range of the middle- 
mass BH in M82, assuming its host corresponds to the 2.2 /v,m secondary peak cluster. The upper 
limit of the middle-mass BH mass and the 2.2 jim secondary peak cluster mass have been estimated 
from IMF calculations, and the lower limit of the BH mass has been estimated from the hard X-ray 
observations. Cross marks indicate the data taken from Magorrian et al. (1998), and a solid line is 
a linear fitting of their data. We did not include their upper limit data in this diagram. This figure 
clearly shows that the data of M82 is on the trend of the mass of the supermassive black holes and 
those of the host bulges. 



This figure "figurel.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/00 11 07 1 v 1 



This figure "figure2.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/00 11 07 1 v 1 



This figure "figure3.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/00 11 07 1 v 1 



